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Abstract 

Ab initio molecular dynamics simulations are used to investigate the structure 
and electronic properties of the liquid Ag-Se system at three compositions. 
The realism of the simulations is demonstrated by comparison with diffrac- 
tion data for the stoichiometric case Ag2Se. As the Se content is increased 
beyond the stoichiometric value, short-lived Se n complexes are formed. The 
concentration of complexes and the associated changes of electronic structure 

can be explained using a simple ionic model. 
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Many binary liquids show a dramatic dependence of their structure and electronic prop- 
erties on composition [|1|]. A celebrated example is the Cs-Au system, in which the pure 
elements are excellent metals, but nevertheless the equiatomic mixture has a very low con- 
ductivity and a structure characteristic of a molten salt 0. These effects arise from the 
electronegativity difference between the elements, and the resulting charge transfer, partial 
ionicity and atomic ordering. In systems where one of the elements is a semiconductor in 
the liquid state - for example alloys of metals with S or Se - even richer behavior can be 
expected, since variation of composition should change the bonding from metallic through 
partially ionic to covalent. We report here a set of simulations of the Ag-Se system per- 
formed using ab initio molecular dynamics (AIMD), which we have used to explore these 
effects. 

The Ag-Se system has been studied by a variety of experimental techniques [|[, but 
understanding of its structure is limited, since diffraction measurements have been made 
only near the stoichiometric composition Ag 2 Se |3J], and these yield only the total structure 
factor and not the partial correlations between the two elements. The present simulations 
allow us to build up a detailed picture of the composition-dependent structure, and also 
give important new insights into the electronic structure and into dynamical aspects of 
the bonding, which would be hard to probe experimentally. A number of previous AIMD 
simulations of liquid metals and semiconductors have been reported || , but the Ag-Se system 
presents a significant technical challenge, as will become clear. 

We shall show that our simulations are in excellent agreement with available structural 
measurements at the stoichiometric composition Ag 2 Se, where we find that a partially ionic 
description is appropriate. As soon as the Se content is increased beyond the stoichiometric 
value, a major structural change begins, consisting of formation of Se n complexes. This 
is accompanied by changes of the electronic density of states (DOS) associated with the 
formation of covalent Se-Se bonds. Both the concentration of Se n complexes and the changes 
in the DOS can be understood using a simple ionic model. In spite of their major structural 
role, the complexes turns out to be surprisingly short-lived, the typical bond lifetime being 
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~ 0.5 ps. 

Our AIMD technique, like the methods pioneered by Car and Parrinello ||, uses den- 
sity functional theory, pseudopotentials, and a plane- wave representation of the orbitals 
However, instead of treating the electronic degrees of freedom as fictitious dynamical vari- 
ables, we use conjugate-gradient minimization || to reach the self-consistent ground state at 
each step. The Hellmann-Feynman forces are then used to integrate the classical equation of 
motion of the ions. To handle the semi-metallic nature of the system, we use Fermi-surface 
smearing, with the electronic occupation numbers treated as auxiliary dynamical variables 
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The details of our calculations are as follows. We use ab-initio norm-conserving non-local 
pseudopotentials, with exchange and correlation included via the local density approxima- 



tion in Ceperley- Alder form [12]. The Ag pseudopotential has been optimized using the 



method of Lin et al. ||13|| , which is a refinement of the scheme due to Rappe et al. [14 



The Se pseudopotential does not require optimization and the standard Kerker method JTo 
suffices. The s and p components of the Ag pseudopotential were generated using the atomic 
configuration 4c? 10 5s ' 25 5p 0,25 , and the d component using the configuration 4d 10 5s 5 . The 
core radii were 2.0, 2.0 and 2.5 a.u. for the s, p and d components respectively. For Se we 
used 4s 2 4p 4 for the s and p waves and 4s 2 Ap 275 3d 0,25 for the d wave; the core radii for 
were chosen to be 2.0, 2.0, and 2.3 a.u. for the s, p and d components respectively. We 
use the pseudopotentials in Kleinman-Bylander separable form |16|] with the s-wave treated 



as local; the non-local parts of the pseudopotentials are treated in real space fl7fl . Tests 



on the low-temperature crystal phase of Ag 2 Se show that these pseudopotentials reproduce 
the experimental structure very well and that a plane-wave cutoff of 400 eV is sufficient 
to achieve good convergence of the total energy; this cutoff is used throughout the present 
work. We have also tested these potentials on AgCl flI8| , pure Se and GaSe. In every case 



we find very good agreement with the experimental structure. 

AIMD calculations on systems containing post-transition metals such as Ag are extremely 
demanding, because it is essential to include the (^-electrons as valence electrons, and because 



of the large plane- wave basis set needed to represent the <i-orbitals. As a result, very few such 
AIMD simulations have been reported Jl9|,^0j. The present work relies heavily on parallel 



computation methods, and has been performed using the parallel ah initio code CETEP |2T| 
running on a Cray T3D machine. 

Our simulations of the Ag-Se liquid alloys have been performed on a system of 69 atoms 
in a cubic box with the usual periodic boundary conditions. We use T-point sampling, 
and a Fermi-smearing energy width of 0.2 eV. The Verlet algorithm is used to integrate 
the ionic equation of motion, with a time step of 3 fs. We have performed simulations at 
the temperature T ~ 1350 K for three concentrations of Agi^Se^, namely: x = 0.33 (46 
Ag atoms and 23 Se atoms), x = 0.42 (40 Ag atoms and 29 Se atoms) and x = 0.65 (24 
Ag atoms and 45 Se atoms). The simulations are performed at a density which is linearly 
interpolated between experimental values for £-Ag 2 Se and £-Se. 

To initiate the simulations, we exploit the fact that an empirical pair-potential model 
|22|| has been developed for the stoichiometric Ag 2 Se system, which gives a reasonable re- 
production of the liquid structure when used in classical molecular dynamics. We begin by 
making simulations with this empirical model, and we then switch over to AIMD and let 
the system equilibrate for a further 1 ps before collecting data over the next 2 ps. We reach 
the other Agi^Se^ compositions by replacing some of the Ag atoms by Se atoms and then 
equilibrating for 1 ps at the new composition; production runs of 2 ps are again performed 
in each case. 

Fig. [I] shows the neutron-weighted static structure factor S(k) obtained from our simu- 
lation for the Ag 2 Se composition, compared with experimental data at T = 1150 K [Q. The 
overall agreement between theory and experiment is excellent, with all the main features 
having the correct position and magnitude. It is also remarkable that our simulation repro- 
duces the position of the small pre-peak at wave vector k ~ 1.7 A ; its height, however, 
is lower than the experimental value. This discrepancy is probably due to the finite size of 
our simulation cell and the resulting lack of fc-space resolution. Our simulation allows us to 
analyze the origin of the peaks in S(k) in terms of the calculated partial structure factors 



SAg-Ag(^), SAgSe(k) and Sse-Se(k). We find that the main peak in S(k) at k ~ 2.7 A -1 is due 
to peaks in S , A g -A g (^) and £>Ag-Se(fc), which reinforce each other. The pre-peak in S(k) arises 
from a positive peak in Ss e -Se(k) and a negative peak in S\ g _s e {k), which partially cancel 
each other. A more detailed analysis of the structure factors will be given elsewhere. 

The partial pair correlation functions g\ g -A g {r), g\ g -Se(r) and gs e -Se(r) for the three con- 
centrations are displayed in Fig. |]. The results show that increase of Se content causes 
dramatic changes in gse-Se(r)- As Se content increases, a short- distance peak builds up. At 
stoichiometry, gs e -Se( r ) exhibits broad peaks at 3.99 and 4.83A, and only a very weak tail 
below 3.0 A. For x = 0.42, instead of a tail, there is a short- distance peak at 2.35 A. The 
position of the main peak has shifted to 4.72 A and has decreased in magnitude. The peak 
around 4.0 A seen at stoichiometry has merged with the main peak to give rise to a shoul- 
der. At the last concentration, the short- distance peak is dominant and gse-Se{f) shows little 
structure beyond 4.0 A with a low and broad second peak. We note that the radius 2.35 A 
associated with the short- distance peak is very close to the Se-Se covalent bond length in 



crystalline and liquid Se ||23|| . The growth of the short- distance peak can be characterized by 
the Se-Se coordination number within a sphere of radius 3.0 A, which we calculate to be 0.1, 
0.71 and 1.7 for x = 0.33, 0.42 and 0.65 respectively. By contrast with the major changes 
in f^se-Se^), the other two correlation functions change rather little with composition. 

More insight into the structural changes can be gained from 'snapshots' of the ionic 
positions at the three concentrations (Fig. [3]). To aid the eye, bonds have been drawn 
between Se atoms separated by less than 3.0 A. As expected from the results for g a p(r), the 
major feature is the formation of Se clusters as x exceeds 0.33. Already at x = 0.42, the Se 
atoms bond not only into dimers but also into larger Se n clusters. At this composition, only 
~ 48 % of the Se atoms are bonded, 76 % of the bonded atoms being one-fold coordinated 
and 22 % two-fold, the remaining 2 % having higher coordination. The dominance of one- 
fold and two-fold coordination means that most of the clusters are either dimers or Se n 
chains. 
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At the composition x = 0.65, most of the Se atoms are in clusters, with only 7 % being 
isolated. The proportions of bonded Se in one-fold and two-fold coordination are now 35 % 
and 40 %, with a non- negligible 15 % being three-fold coordinated. The liquid is thus 
composed of one-dimensional Se ra chains interconnected via three-fold coordinated atoms, as 
shown in Fig. |3]c. The Se n clusters at this composition are large, with n typically > 20. 

The structural changes are intimately linked to the electronic structure. The calculated 
electronic density of states (DOS) and the local DOS on Ag and Se atoms for the three 
compositions are shown in Fig. 4. The main features are Se(4s) states at -12 eV, Ag(4<i) 
states at -4 eV and Se(4p) states in the region above -7 eV; Ag(5s-p) states extend upwards 
from roughly the Fermi level E-p, and there is a significant hybridization between these 
and the Se(4p) states. The Se(4p) band of states stands out clearly from other parts of 
the DOS at all three compositions. For Ag 2 Se, E F lies in a pseudo-gap at the top of this 
band; the number of occupied 4p states per Se atom is equal to 3, and a partially ionic 
model approximating to Ag 2 f Se = is appropriate. As x increases beyond 0.33, the number 
of occupied 4p states per Se atom falls below 3, and we attribute this to the formation 
of unoccupied anti-bonding 4p combinations associated with clusters. These anti-bonding 
states are hybridized with Ag(5s-p) states, and only become clearly visible at x = 0.65. 
We also note that with increasing x the Se(4s) band broadens strongly, and we can show 
that this arises from the formation of bonding and anti-bonding Se(4s) combinations in the 
clusters. 

Our analysis of the structure and the DOS thus leads to a simple approximate picture of 
liquid Ag-Se alloys. At the Ag 2 Se composition, the system consists approximately of isolated 
Ag + and Se = ions. Our picture is that as the Se content increases, chain-like Se n clusters are 
formed, and the system consists of a mixture of Ag + , Se = ions and (Se n ) = complexes. As 
confirmation of this picture, we can use charge-balance arguments to estimate the number 
of Se-Se bonds at any composition. If the system consisted entirely of Ag + and Se = ions, 
the net charge on a system of iV Ag and N Se ions would be N^ g — 2N Se . If the clusters 
are all (Se n ) = , the formation of every Se-Se bond reduces the charge by two units. For 
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electroneutrality, the number of bonds must therefore be N$ e — \N^ g and the coordination 
number of the short-distance peak in gs e -Se(f) must be 2 — NA g /Ns e = 3 — 1/x. This 
gives predicted values of 0, 0.62 and 1.46 for the selenium-selenium coordination number at 
the three compositions, which are quite close to the values 0.1, 0.71 and 1.7 reported above. 
There are experimental indications for the occurrence of complex ions in other chalcogenide- 



metal alloys, for example (Se2) pairs in CuSe [g4J] and (Te n ) chains in the Te rich side of 
K-Te |5|. 

In view of the dominant structural role played by Se n clusters, the dynamics of their 
formation and dissolution is of great interest. We have performed a statistical analysis on 
the lifetime of Se-Se bonds, which shows that under the conditions we have studied they 
have a very short half-life on the order of 0.5 ps. A striking effect that seems to be linked 
with this is that the Se diffusion coefficient increases with increasing Se content. Details of 
these dynamical effects will be described elsewhere. 

In conclusion, we have performed AIMD simulations on the liquid Ag-Se system to in- 
vestigate how the atomic ordering and electronic structure evolve with composition in this 
typical alloy between metallic and semiconducting elements. The close agreement with 
neutron-diffraction data confirms the realism of the simulations. We have shown that the 
formation of chain-like Se n complexes begins as soon as the Se content exceeds the stoichio- 
metric value, and that the covalent Se-Se bond-length is very close to its value in crystalline 
and liquid Se. The concentration of Se-Se bonds and the changes of electronic structure can 
be understood using a simple ionic model in which the liquid is a mixture of Ag + , Se = and 
(Se n ) = ions. Although the Se„ complexes are structurally so important, individual Se-Se 
bonds are extremely transient, with a lifetime of less than a picosecond. 

ACKNOWLEDGMENTS 

The work was done within the U.K. Car-Parrinello Consortium, which is supported by 
the High Performance Computing Initiative. A time allocation on the Cray T3D at EPCC 



7 



is acknowledged. Discussions with J. Enderby and A. Barnes, and technical help from I. 
Bush, M. Payne, A. Simpson and J. White played a key role. JMH's work is supported by 
EPSRC grant GR/H67935. The work used distributed hardware provided by EPSRC grants 
GR/H31783 and GR/J36266. 



8 



REFERENCES 

[1] J. E. Enderby and A. C. Barnes, Rep. Prog. Phys. 53, 85 (1990). 
[2] F. Hensel, Adv. Phys. 28, 555 (1979). 

[3] H. Endo, M. Yao, and K. Ishida, J. Phys. Soc. Japan 48, 235 (1980); V. M. Glasov, S. 
M. Memedov, and A. S. Burkhanov, Sov. Phys.-Semicond. 20, 263 (1986); S. Ohno, A. 
C. Barnes, and J. E. Enderby, J. Phys.: Condens. Matter 6, 5335 (1994). 

[4] D. L. Price, M.-L. Saboungi, S. Susman, K. J. Volin, J. E. Enderby, and A. C. Barnes, 
J. Phys.: Condens. Matter 5, 3087 (1993). 

[5] I. Stich, R. Car, and M. Parrinello, Phys. Rev. Lett. 63, 2240 (1989); G. Galli and M. 
Parrinello, J. Chem. Phys. 95, 7504 (1991); G. A. de Wijs, G. Pastore, A. Selloni, and 
W. van der Lugt, Europhys. Lett. 27, 667 (1994); M. Schone, R. Kaschner, and G. 
Seifert, J. Phys.: Condens. Matter 7, L19 (1995). 

[6] R. Car and M. Parrinello, Phys. Rev. Lett. 63, 2471 (1985). 

[7] See e.g. G. P. Srivastava and D. Weaire, Adv. Phys. 36, 463 (1987); J. Ihm, Rep. Prog. 
Phys. 51, 105 (1988); M. J. Gillan, in Computer Simulation in Materials Science, eds. 
M. Meyer and V. Pontikis, p. 257 (Kluwer, Dordrecht, 1991); G. Galli and M. Parinello, 
ibid, p. 283. 

[8] M. C. Payne, M. P. Teter, D. C. Allan, T. A. Arias, and J. D. Joannopoulos, Rev. Mod. 
Phys. 64, 1045 (1992). 

[9] M. J. Gillan, J. Phys.: Condens. Matter 1, 689 (1989). 

[10] G. Kresse and J. Hafner, Phys. Rev. B 49, 14251 (1994). 

[11] M. P. Grumbach, D. Hohl, R. M. Martin, and R. Car, J. Phys.: Condens. Matter 6, 
1999 (1994). 

[12] J. Perdew and A. Zunger, Phys. Rev. B 23, 5048 (1981). 

9 



[13] J.-S. Lin, A. Qteish, M. C. Payne, and V. Heine, Phys. Rev. B 47, 4174 (1993). 

[14] A. M. Rappe, K. M. Rabe, E. Kaxiras, and J. D. Joannopoulos, Phys. Rev. B 41, 1227 
(1990). 

[15] G. P. Kerker, J. Phys. C 13, L189 (1980). 

[16] L. Kleinman and D. M. Bylander, Phys. Rev. Lett. 48, 1425 (1982). 

[17] R. D. King-Smith, M. C. Payne, and J. S. Lin, Phys. Rev. B 44, 13063 (1991). 

[18] F. Kirchhoff, J. M. Holender, and M. J. Gillan, Phys. Rev. B 49, 17 420 (1994). 

[19] A. Pasquarello, K. Laasonen, R. Car, C. Lee, and D. Vanderbilt, Phys. Rev. Lett. 69, 
1982 (1992). 

[20] G. Kresse and J. Hafner, Phys. Rev. B 48, 13115 (1993). 

[21] L. J. Clarke, I. Stich, and M. C. Payne, Comp. Phys. Comm. 72, 14 (1992). 

[22] J. P. Rino, Y. M. M. Hornos, G. A. Antonio, I. Ebbsjo, R. K. Kalia, and P. Vashishta, 
J. Chem. Phys. 89, 7542 (1988). 

[23] K. Suzuki, Ber. Bunsenges. Phys. Chem. 80, 689 (1976). 

[24] A. C. Barnes and J. E. Enderby, Phil. Mag. B 58, 497 (1988). 

[25] J. Fortner, M.-L. Saboungi, and J. E. Enderby, Phys. Rev. Lett. 69, 1415 (1992). 



10 



FIGURES 

FIG. 1. The total neutron weighted structure S(k) factor of Ag2Se. Solid line and circles 
represent simulation and experimental results Q respectively. 

FIG. 2. Partial radial distribution functions g a p{r) of £-Agi_ x Se x at concentrations 2=0.33 
(full line), 2=0.42 (dotted line) and 2=0.65 (dot-dashed line). 

FIG. 3. Snapshots of typical configurations of £-Agi__ 1 .Se 1 . at concentrations (a) x=0.33, (b) 
2=0.42 and (c) 2=0.65. Silver atoms are shown as black spheres, selenium atoms as gray spheres. 
Bonds are drawn between Se atoms with separation < 3.0 A. Bonds to atoms in neighboring cells 
are represented with two-colored sticks. 

FIG. 4. Density of states (solid curve) and local densities of states (LDOS) for Ag (chain curve) 
and Se (dotted curve) from simulations of Agi— ^Se^ at 2 = 0.33, 0.42 and 0.65. For clarity, the 
scale used for the Se LDOS is four times that used for the Ag LDOS. The vertical dotted line 
marks the Fermi energy. 
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